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[i] The low latitude boundary of the proton aurora (known as the Isotropy Boundary 
or IB) marks an important boundary between empty and full downgoing loss cones. 

There is significant evidence that the IB maps to a region in the magnetosphere where the 
ion gyroradius becomes comparable to the local field line curvature. However, the location 
of the IB in the magnetosphere remains in question. In this paper, we show simulated 
proton precipitation derived from the Field Line Curvature (FLC) model of proton 
scattering and a global magnetohydrodynamic simulation during two substorms. 

The simulated proton precipitation drifts equatorward during the growth phase, 
intensifies at onset and reproduces the azimuthal splitting published in previous studies. 

In the simulation, the pre-onset IB maps to 7-8 R E for the substorms presented and the 
azimuthal splitting is caused by the development of the substorm current wedge. 

The simulation also demonstrates that the central plasma sheet temperature can 
significantly influence when and where the azimuthal splitting takes place. 

Citation: Gilson, M. L., J. Raeder, E. Donovan, Y. S. Ge, and L. Kepko (2012), Global simulation of proton precipitation due 
to field line curvature during substorms, J. Geophys. Res., 117, A05216, doi:10.1029/2012JA017562. 


1. Introduction 

[ 2 ] Much effort has been spent to determine the Earth’s 
magnetic field topology under various inputs from the solar 
wind. This infonnation is useful because it can be used to 
draw connections between auroral features and magneto- 
spheric phenomena. In a recent example of this, Nishimura 
et al. [2010] was able to correlate individual patches of 
patchy aurora with chorus waves in the magnetosphere. In 
other recent papers, authors have speculated about the origin 
of the wavy structure that appears along the substonn onset 
arc a few seconds before onset [e.g., Raeder et al., 2010], 
For substorms, much debate has existed about where the 
breakup arc originates compared to the location of the near- 
earth reconnection. Since the mapping is uncertain, this 
remains a fundamental disagreement between researchers. 

[ 3 ] One of the difficulties associated with mapping auroral 
features during substorms arises because the magnetotail is 
highly stretched and can globally reconfigure over short peri- 
ods of time. Because typical empirical models [ Tsyganenko , 
1989, 2002a, 2002b] are parameterized by geomagnetic indi- 
ces and Interplanetary Magnetic Field (IMF) measurements 
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(K P , DST, AE, B z , etc.) they only represent an average field 
configuration and are therefore most applicable to general 
statistical studies as opposed to particular events. To circum- 
vent this problem, authors often use data from low-altitude 
satellites in the auroral oval to leam about particle distributions 
and field structure in the source regions [e.g., Sergeev et al., 
1993,2012], 

[ 4 ] One of the more promising routine measurements 
made by low altitude satellites is the location of the isotropic 
boundary (IB). Equatorward (earthward) of the low latitude 
boundary of the proton aurora, low altitude satellites mea- 
sure double loss cone ion distributions indicating bounce- 
trapped particles. However, in measurements at altitudes 
above the proton auroral regions, satellites typically measure 
a full downgoing loss cone [Donovan et al., 2003a], The IB 
is the location of the transition between the full and empty 
downgoing loss cones. The particles in the downgoing loss 
cone have energies consistent with central plasma sheet 
(CPS) ions [Eather, 1967] and it is now commonly accepted 
that they are generated by scattering in the CPS. Generally, 
the IB is a function of particle energy and species, however, 
it generally corresponds well with b2i boundary of Newell 
et al. [1996], defined as the location of the maximum of 
the integrated ion energy flux. 

[ 5 ] The IB and the b2i have a long history of being used to 
infer the degree of magnetic field stretching in the inner 
magnetotail. Simulations have shown that in regions where 
the field line curvature (FLC) becomes comparable to the 
particle gyroradius there is significant pitch angle scattering 
[Sergeev and Tsyganenko, 1982], The criterion for strong 
enough scattering to fill the loss cone found by Sergeev et al. 
[1983] is that the magnetic field radius of curvature divided 
by the particle gyroradius is less than 8. It is now more 
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common to see the K-parameter [Buchner and Zelenyi, 
1987]. In this form, the strong scattering criterion is: 



where R c is the magnetic field radius of curvature and p is 
the particle gyroradius. 

[6] If the scattering is caused by FLC, then the IB maps to 
an important transition from moderately stretched to highly 
stretched field lines in the magnetosphere. This inspired 
Sergeev et al. [1993] and Newell et al. [1998] to show that 
the magnetic field inclination at geosynchronous orbit is 
highly correlated with the IB and b2i respectively on the 
midnight meridian. They used the auroral boundaries to add 
an additional constraints on the existing empirical mapping 
models. The correlation was good enough (r ~ 0.9) to lead 
Gvozdevsky and Sergeev [1996] to define a magnetotail 
(MT) index based on the IB latitude. More recently, Meurant 
et al. [2007] showed that the IB remains well correlated with 
the extent of magnetic field stretching even during active 
times like the substorm expansion phase. They also showed 
that the stretching typically decreases first in the dawn sector. 

[ 7 ] In addition to correlating with the degree of stretching 
in the magnetotail, the onset arc typically sits near the pole- 
ward boundary of the proton precipitation [Lessard et al., 
2007; Deehr and Lummerzheim, 2001]. Thus, the onset arc 
definitely maps tailward of the IB during substorms. As 
much debate revolves around the timing of the auroral 
brightening compared to the transit time of observed sub- 
storm signatures in the tail, the location of the IB could at 
least place a lower bound on the mapping of the onset loca- 
tion. Synoptic descriptions of how the electron and proton 
auroral precipitation boundaries are positioned relative to 
one another as a function of magnetic local time have been 
presented in Creutzberg et al. [1988] and Mende [2003]. 

[ 8 ] Because the importance of the IB and the success of 
the FLC model, observations of the IB and b2i have also 
been inferred from meridian scanning photometers (MSP) 
[Donovan et al., 2003b] and from global images taken by the 
SI- 12 instrument on the IMAGE spacecraft [Blockx et al., 
2005; Meurant et al., 2007]. However, despite the large 
number of observations, where the IB maps in the tail is still 
an open question. 

[ 9 ] Recently, Gilson et al. [2011] focused further on the 
longitudinal evolution the proton aurora during substorms 
and showed that the proton aurora commonly splits after the 
onset with eastward and westward traveling precipitation 
fronts. The splitting is more likely to occur for substorms 
with lower AL. Based on the FLC model, they speculated 
that the split region corresponds to the edges of the substorm 
current wedge (SCW) as the field curvature inside the SCW 
should be much larger than outside. However, in that study, 
it was unclear why the proton aurora should split during 
some events and not others. One hypothesis was that the 
scattering could continue in regions with a quasi-dipolarized 
magnetic field geometry if the plasma temperature was signif- 
icantly high. This hypothesis was supported by the observation 
that weaker substorms are less likely to split than their stronger 
counterparts because the weaker substorms are probably 
associated with less intense dipolarization. Alternatively, 


wave-particle interactions [Ashour-Abdalla and Thorne, 
1978] could be responsible for scattering particles into the 
loss cone in the SCW. 

[ 10 ] A study by Donovan et al. [2003a] showed that the 
energy dependence of the IB latitude is rarely completely 
consistent with what is expected from the FLC model. It is 
reasonable to expect that higher energy particles will scatter 
in a less stretched field due to their larger average gyroradii. 
This implies that the IB for high energy particles should be 
observed at lower latitudes than the IB of low energy parti- 
cles to be consistent with the FLC model. Of approximately 
1000 FAST passes through the auroral oval, only ~ 1 0% of 
the transits had energy orderings completely consistent with 
the FLC model while ~10% were completely inconsistent. 
The discrepancy could be due to wave-particle interactions. 
If that turns out to be the case, using the proton aurora to 
mark the low k boundary may not valid. However, it is also 
possible that some other physical process is occurring at the 
inner edge of the CPS to reverse the ordering. One candidate 
is charge exchange leading to precipitating heavier ions 
[Kistler et al., 1998] since the FAST -ESA instrument used in 
the Donovan et al. [2003a] study did not have the mass 
resolution to separate out that effect. 

[ 11 ] In this study, we use the FLC model with the 
OpenGGCM global MHD code to simulate the proton pre- 
cipitation during substorm conditions using realistic solar 
wind input. We will present two substorms. The first is the 
March 23, 2007 “THEMIS First Light” substorm analyzed 
by Angelopoulos et al. [2008], Raeder et al. [2008], Runov 
et al. [2008], Keiling et al. [2008], Lessard et al. [2009] 
and others. The second is a substorm observed by the 
IMAGE Far Ultraviolet (FUV) camera [Mende et al., 2000b, 
2000a, 2000c] on April 28, 2001. In the simulations, the 
proton precipitation splits longitudinally and forms west- 
ward and eastward traveling fronts similar to the description 
given by Gilson et al. [201 1], This splitting is well correlated 
with the longitudinal growth of the SCW in the simulation. 
Similar splitting is seen by the IMAGE FUV camera for the 
April 28, 2001 substorm. Therefore, we show that the FLC 
model is sufficient to explain the observed proton auroral 
splitting. We also show that for some conditions the pre- 
cipitation due to FLC can exist in a quasi-dipolarized SCW 
which eliminates the need to consider wave particle inter- 
actions to explain the lack of proton auroral splitting during 
weaker substorms. 

2. OpenGGCM 

[ 12 ] To address the mapping problem, we are using 
OpenGGCM for this study. OpenGGCM is a global simu- 
lation which solves the MHD equations on a stretched Car- 
tesian grid. The advantage of this grid structure is that the 
resolution is able to be increased in areas where it is needed 
(e.g., the magnetotail) while keeping the total number of 
cells manageable. The main advantage of using a global 
model for this study is that auroral features can be mapped to 
structures in the tail self-consistently without needing to wait 
for a fortunate conjunction between satellites and/or ground 
observations on a clear day. Also, in some cases, the simu- 
lation mapping can explain observed data better than avail- 
able empirical models [Angelopoulos et al, 2008]. 
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Figure 1. Synthetic keogram along the Gillam meridian for 
the March 23, 2007 simulation, (top) The electron precipita- 
tion accelerated in regions of upward field aligned current, 
(middle) electron precipitation calculated from the thermal 
electron flux, and (bottom) the proton precipitation calcu- 
lated as described in the text. The equatorward drifting and 
intensity increase near the onset is evident. 


[ 13 ] The boundary conditions on the dayside face of the 
simulation domain are set by data from a solar wind monitor 
and the other five faces have open boundary conditions. The 
three dimensional structure of the solar wind data is inferred 
by the procedure described in Raeder et al. [2001, section 3] 
using the minimum variance analysis of Sonnerup and 
Cahill [1967]. Hereafter, that procedure is referred to as 
MINVAR. At times, the solar wind may not be ordered in 
simple sheets and so the MINVAR procedure may not 
accurately reproduce the field. In those cases, if the B x 
component is not varying much, setting B x to the average 
value of the interval can be a reasonable approach. In other 
situations, more complicated methods may be used to obtain 
the dayside boundary condition [e.g., Russell et al., 2000]. 

[ 14 ] The inner boundary is calculated by mapping simu- 
lation current and electric field to the ionosphere and using 
empirical models to calculate the electron precipitation 
parameters and potential drop. The simulation calculates the 
thermal electron flux from the inner magnetosphere due to 
scattering and also energy flux from electrons accelerated in 
regions of upward field aligned current. The second popu- 
lation likely corresponds to what is commonly referred to as 
discrete aurora whereas the first corresponds to diffuse 
aurora. The CTIM (Coupled Thermosphere-Ionosphere 
Model) module uses these parameters to compute the iono- 
spheric conductances which are then used to solve for the 
ionospheric potential. See Fuller-Row ell et al. [1996] for 
more the details of CTIM and Raeder et al. [2001] for details 
on the coupling between OpenGGCM and CTIM. For more 
information on OpenGGCM, see Raeder [2003] and Raeder 
et al. [2008]. 

3. Method 

[ 15 ] To calculate the precipitating proton energy flux in 
OpenGGCM, field lines are integrated using a Runge-Kutta 


algorithm with foot points on the ionospheric grid. The 
typical resolution for the ionospheric grid is half a degree in 
latitude and three degrees in longitude. If the field line is 
closed, the average particle gyroradius (p) and field radius of 
curvature ( R c ) (and therefore n) are calculated in the CPS 
from the MHD variables assuming a simple Maxwellian 
particle distribution. In other words, 

s/2mik b Ti 

<2) 

and the radius of curvature is calculated by: 

~=-(b-V)b (3) 

where 7] and m L are the ion temperature and mass, k h is the 
Boltzmann constant, e is the fundamental unit of charge, B is 
the magnetic field and b is a unit vector in the direction of B. 
Since the loss cone is typically full and k is known, the 
diffuse precipitating energy flux is calculated using 

F =f(n)n i k h T i sj 2k b Ti/irmi (4) 

where », is the ion number density. Equation (4) is the same 
used by Raeder et al. [2008] for diffuse electron precipita- 
tion with parameters for ions substituted and with an addi- 
tional dependence on k. For this paper, /( k) was taken to be 
a simple rectangle function: 

/(k) = / 1 : 0 < /c < x/8 (5) 

( 0 : otherwise 

The simple form of f[n) in equation (5) was chosen because 
of the sharp cutoff in pitch angle diffusion reported by 
Sergeev et al. [1983]. Other forms of J[k) have been reported 
in the literature. For example, Delcourt and Martin [1994] 
and Delcourt et al. [1996] showed that in the regime where 
k < 1 , the scattering does not efficiently move particles into 
the loss cone. Liu et al. [2007] used /(k) 00 B z exp(— 0.97 v / k) 
to explain the pre-onset fading of the proton aurora. However, 
it is unclear how well the MHD approximation resolves the 
field curvature (and therefore k) in the highly stretched CPS 
as one of the assumptions of MHD is that the average gyro- 
radius is small compared to the average length scale. There- 
fore, these forms have not been used for this study. 

[ 16 ] For each substorm presented here, the proton precip- 
itation is poleward of the discrete electron precipitation 
duskward of the onset and equatorward of the discrete 
electron precipitation on the dawnward side consistent with 
the statistical studies by Creutzberg et al. [1988] and Mende 
[2003]. The proton precipitation also intensifies and expands 
poleward in the simulation. 

4. Event 1: March 23, 2007 

[ 17 ] March 23, 2007 was one of the first substorms 
observed by the THEMIS probes, and therefore the aspects 
of this event have been studied in detail by Angelopoulos 
et al. [2008], Raeder et al. [2008] and many others. The 
latter simulated the event using the OpenGGCM code. For 
specific details of that simulation and analysis of some of the 
results, the reader is referred to Raeder et al. [2008]. The 
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Figure 2. (a-d) Simulated proton aurora mapped along field lines to the furthest point with contours of 
B z at different times as marked, (e-h) B z with contours provided for reference. B z and contours extracted 
from the plane of maximum plasma (3 as a proxy for the central plasma sheet. 


simulation mapping agreed better with the observations 
than the TSOI mapping [Tsyganenko, 2002a, 2002b] 
because the TSOI mapped locations of the THEMIS probes 
lie outside the westernmost location of the auroral activation 
even though all the probes recorded substorm signatures 


(dipolarization, flows and energetic particle flux increases). 
A synthetic keogram of the simulated aurorae is presented in 
Figure 1. From the keogram, it is obvious that the simulation 
captures the equatorward drifting during the growth phase 
commonly associated with substorms. Plots from four 
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Figure 3. Simulated proton aurora (color) and discrete electron precipitation (contours) at different times 
in the simulation of the March 23, 2007 substorm. Local noon is at the top of each plot. The proton pre- 
cipitation splits on the eastward edge of the developing westward traveling surge. 


different times in the early stages of the substorm are shown 
in Figure 2. The first column shows the simulated proton 
precipitation mapped along magnetic field lines to the CPS. 
For the purposes of data visualization, the magnetic field is 
also plotted on the CPS (defined by the surface of maximum 
plasma (3). The simulation shows a Dipolarization Front 
(DF) propagating earthward from the reconnection site at 
about —15 R e prior to the auroral brightening. Once the 
front reaches the inner magnetosphere, the flux piles up and 
grows toward the flanks forming the SCW (Figure 2). 

[is] The proton precipitation shows the signature of the 
SCW in the ionosphere. Since the FLC model requires a 
highly stretched (tail-like) magnetic field or large gyroorbits 
(high temperature) to scatter particles into the loss cone, the 
dipolarized region can be seen in sudden latitudinal increa- 
ses in the IB location causing a longitudinal split in the 


proton precipitation (compare Figures 2a-2d with Figure 3). 
Another possibility is that the split is caused by the longi- 
tudinal dependence of the temperature. This would require a 
region of low temperature plasma at or near the onset 
meridian with high temperature plasma on either side. The 
simulation of this event shows the opposite. High tempera- 
ture plasma is injected from the reconnection site into the 
dipolarized region and expands outward with (but inside) the 
current wedge. 

[ 19 ] The simulation agrees with the interpretation of the 
longitudinal splitting of Gilson et al. [2011], however, there 
are unfortunately very limited ground or space proton aurora 
data for this event to compare with the simulation. Also, this 
simulation does not explain why the proton aurora splits in 
stronger events but not in weaker events. Because of the 
limited data, we use the IMAGE data set to find an event 
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Figure 4. Provisional AE and DST indices for April 28, 2001. The vertical black line is the approximate 
breakup time as seen by the IMAGE spacecraft. 


which clearly shows the splitting to gain more insight into 
what causes it. 

5. Event 2: April 28, 2001 

[ 20 ] On April 28, 2001, IMAGE recorded a substorm with 
an initial brightening shortly after 13:00 UT. IMAGE started 
collecting data on a pass over the northern hemisphere at 
about 11:05 UT. The first few images show significant 
westward propagating aurora, probably due to a previous 
substorm. After 11:54 UT, the images show almost no 
activity other than the dimly illuminated auroral oval until 
the breakup arc begins to brighten at the onset location at 
13:07:49 UT [Frey et al, 2004], This brightening is 
accompanied by a significant enhancement in the AL index, 
poleward expansion and westward traveling motion of the 
electron precipitation indicating that this is a true substorm 
and not a pseudo-breakup. The AU, and AL indices are 
plotted in Figure 4. For this event, the proton aurora also 
brightens significantly and expands poleward. 

[ 21 ] The magnetosphere was in a reasonably disturbed 
state with K P between four and six and DST ranging from 
40 nT to about —35 nT starting at 06:00 UT until about 
17:00 UT (Figure 4). However, at the time of the substorm 
onset the DST was around zero. The energy was loaded into 
the magnetosphere by an above average solar wind speed 
(V x 650 km/s) and an average number density (n ~ 5 cc -3 ). 
At the onset of the substorm, a sharp increase in the dynamic 
pressure was brought on by an increase in the density to 
approximately 8.5 cc - ’. The y-component of the IMF was 
significantly above average at 1 5 nT throughout the period of 
interest while the z-component varied greatly. Approxi- 
mately one hour prior to the substorm, B z turned southward 
with a magnitude approaching 10 nT for about half an hour. 


It then turned northward with a magnitude of about 1 5 nT and 
proceeded to decrease to about — 5 nT at the time of onset (see 
Figures 5, 6 and 9). 

[ 22 ] Figure 7 shows the early development of the expan- 
sion phase using the images from the W1C (electron aurora) 
and SI- 12 (proton aurora) cameras on the IMAGE spacecraft. 
Two frames after the initial brightening (at 13:11:54 UT), the 
SI- 12 instrument shows the proton precipitation peak split 
just eastward of the initial brightening at about 23 magnetic 
local time (MLT). The two proton precipitation fronts sepa- 
rate until the gap between them encompasses approximately 
two hours MLT and then remain in this configuration until 
after 14:30 UT. For this event, the eastward precipitation 
front expands westward while its eastward boundary remains 
fixed. The westward front continues to expand westward and 
its eastward boundary remains nearly co-located with the 
eastward boundary of the WTS. 

[ 23 ] After the initial brightening of the aurora in this event, 
a second enhancement of the electron precipitation is evident 
at around 13:44 UT. This corresponds to another enhance- 
ment in the proton precipitation as well as a second drop in 
the AL index. Even with the very active auroral oval at that 
time, the proton precipitation remained split. 

5.1. Simulation and Results 

[ 24 ] For this event, the simulation domain was from 
—500 R e to 24 R e GSE in the x-direction and from —64 R E 
to 64 R e in y and z. We used a high resolution grid with 
nearly 37 million cells (616 x 200 x 300). The minimum 
resolution was approximately 0.2 R E in x, 0.25 R E in y and 
0.16 R e in z. The high resolution region spanned the tail to a 
distance of approximately 40 R E . The grid spacing as a 
function of the grid axis is shown in Figure 8. 
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Figure 5. Solar wind input to the simulation propagated from GSE 224, 16, -24 to the sunward boundary 
of the simulation f ~24 R E upstream). Panels show (from top to bottom) magnetic field (nT), plasma flow 
velocity (km/s), number density (particles per cc), and thermal pressure (pPa). The vertical lines are at 
1 1 :54:00 UT and 13:07:49 UT which are the last significant aurora seen by IMAGE SI-12 and the auroral 
onset respectively. 


[ 25 ] Figure 9 shows the magnetic field measured at the 
Geotail spacecraft located in the magnetosheath. The simula- 
tion does a good job propagating the IMF y and z components 
into the magnetosheath, however the x component is not very 
good. This typically means that the IMF was not well ordered 
in sheets as assumed by the MINVAR procedure. This does 
not have a significant influence on the global dynamics in the 
magnetotail since the dayside reconnection rate is only weakly 
affected by the IMF B x as long as |B X |/B <■ 1 because the 
draping of the IMF around the magnetosphere typically redu- 
ces B x before it interacts with the magnetosphere. 

[ 26 ] The simulation reproduces the equatorward drifting 
of both the proton and electron precipitation maxima during 
the growth phase (Figure 6). Similar to the March 23, 2007 
substorm, the pre-onset IB in the simulation maps to around 
7-8 R e . The onset occurs around 13:04 UT in the simulation. 
However, the typical substorm onset signatures occur over a 
much larger timescale in the simulation than in observations. 
Because of this, the chosen onset time for this event is a 
compromise between the enhancement of a simulated AL 
index and the most explosive brightening of the electron 
precipitation that occurs approximately ten minutes later. At 
13:04 UT, there is slight brightening of the electron precipi- 
tation and the initial poleward expansion accompanied by 
westward propagation (Figure 7, column 3). There is also an 
enhancement of the plasma sheet temperature at that time as 
well which can be inferred from the increased proton energy 
flux in Figure 7. Also, the onset longitude is approximately 
correct at 13:04 UT. If the onset was placed later, the onset 
meridian would be displaced from the actual onset location 
by up to 1.5 hours MLT due to the movement of the WTS. 

[ 27 ] At 13:04 UT, there is already a poleward displace- 
ment in the low latitude boundary (i.e., the IB) of the 


simulated proton precipitation (Figure 7). However, the 
scattering still persists at high latitudes resulting in pseudo- 
split proton precipitation (magenta arrows in Figure 7). 
Unlike the observations, the IB in the pseudo-split region 
gradually moves equatorward preventing the pseudo-split- 
ting from becoming fully split as in the observations. This 
discrepancy will be discussed shortly. Presently the proton 
precipitation in the simulation split again. The second split- 
ting can be seen in Figure 10. This second split clearly 
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Figure 6. Simulated keogram along the 21:30 local time 
meridian during the growth phase of the April 28, 2001 
event. The format is the same as Figure 1. The keogram 
reveals two auroral brightenings. The first is near 12:30 UT 
and the second is near 13:15 UT along this meridian. 
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Figure 7. (a-d) Raw counts from the IMAGE WIC camera (discrete electron precipitation is shown left) 
and from the SI- 12 camera (proton precipitation is shown in the middle) showing the evolution of the pro- 
ton splitting during the first fifteen minutes of the April 28, 2001 substorm. Simulated proton precipitation 
with contours of discrete electron precipitation (shown right) during the late growth phase and into the 
early expansion phase in the simulation. The simulation reproduces some of the early splitting (magenta 
arrow), although it is not as dramatic as in the observations. 


divides the proton aurora into two regions and remains until 
the end of the run. Similar to the observations, this splitting 
has a maximum width of approximately two hours MLT. 

[ 28 ] The proton precipitation for this event is much more 
intense than for the March 23, 2007 substorm which is not 
unreasonable due to the much more extreme solar wind 
conditions. Additionally, the simulation produces a series of 
three substorms (the third substorm is the one at 13:04 UT). 
The first substorm happens at approximately 1 1 :00 UT 


which is consistent with the observations of westward 
propagating aurora. However, it is a very weak substorm and 
it develops very slowly (probably due to the temporal 
proximity with the simulation initial conditions). The second 
substorm occurs around 12:15 UT and is most expanded and 
intense around 12:30 UT. This is somewhat consistent with a 
slight enhancement in the AL (and significant AU 
enhancement) at about 12:00 UT in Figure 4. However, 
there is not a significant auroral signature of this 
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enhancement on the nightside in the WIC data. Therefore, 
the simulation seems to have magnified the magnetic dis- 
turbance into a substorm. Because of this, some of the 
energetic plasma remained in the geosynchronous region at 
the start of the 13:04 UT substorm which accounts for the 


high simulated energy flux at the onset. It also explains why 
the split region does not completely form until later in the 
simulation. As the substorm further energized the plasma, it 
quickly reached a point where the ^-scattering criteria was 
satisfied due to the high average gyroradii. This explains the 
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Figure 9. Comparison of simulation results (solid red line) with Geotail data (dashed green line) located 
in the magnetosheath at around GSE —1 1.5, 25.8, 1.2. The black vertical line is the approximate time of 
the first brightening seen by the IMAGE spacecraft. 
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Figure 10. Simulated proton precipitation with contours of discrete electron precipitation for the April 
28, 2001 substorm. Local noon is at the top. The plot shows the dipolarization front at the poleward 
boundary of the first figure near 2 1 MLT. As that flux is added to the substorm current wedge, the proton 
scattering there is arrested resulting in two distinct proton precipitation regions. 


quasi-splitting in the early expansion phase. It is also likely 
that the simulated magnetic field was over-stretched in the 
expansion phase. This is inferred because if the magnetic 
field was sufficiently dipolarized, the high speed flow bursts 
[Keiling et al., 2009; Birn et al., 2004] and pressure gra- 
dients [Lui, 1996, section 5.6, and references therein] asso- 
ciated with the dipolarization should have created a more 
significant auroral brightening. 

[ 29 ] Figure 11 shows the mapped proton precipitation in 
the left column and B z in the right column for several dif- 
ferent times around the second (complete) splitting of the 
proton precipitation in the simulation. The times are the same 
for the left and right columns and the contours of B z are the 
same to allow easier comparison of features. Similar to the 
March 23, 2007 substorm, this substorm also shows that 


the longitudinal splitting of the proton precipitation is due 
to the strong magnetic field with a high radius of curvature in 
the SCW. However, for this event, the proton precipitation 
does not follow the contours of B z the way it did for the 
March 23, 2007 substorm. Because of this, the reason for the 
splitting in this simulation requires further attention. 

[ 30 ] Shortly after 13:20 UT, a new DF was driven into the 
magnetosphere by a strong flow channel (first panels in 
Figure 13 and Figure 1 1). It is likely that a similar flow burst 
produced the second brightening at 13:44 UT in the obser- 
vations. This DF had low enough curvature to prevent pro- 
ton scattering evidenced by the lack of proton precipitation 
which mapped to it. In the ionosphere, a region void of 
proton aurora which mapped to the DF began to extend equa- 
torward from the poleward boundary at about 22 MLT (first 
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Figure 11 . (a-d) Proton aurora mapped along field lines to the equatorial plane. The contours are of B z in 
the plasma sheet, (e-h) B z on the plane of maximum plasma (3 (i.e., the central plasma sheet). 


plot in Figure 10). The DF strengthened the existing SCW as 
they merged preventing the scattering of protons at the equa- 
torward and poleward boundaries. Thus, the splitting is again 
associated with the increase in curvature and total magnetic 
field in the SCW. The scattering, does not follow well defined 
boundaries in B z for this event because of the strong azimuthal 


temperature dependence. Comparison with Figures 8, 9, and 1 1 
reveals that the scattering in the eastern front is facilitated by 
the high temperature plasma in that region in spite of the mar- 
ginally dipolarized magnetic field there. 

[ 31 ] Assuming a Maxwellian distribution, the mean value 
of k 2 at any given location is directly proportional to the 
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radius of curvature and the magnitude of the magnetic field. 
Since dipolarization typically increases both of these quan- 
tities, the simulation always reproduces a split in the highly 
dipolarized regions of the SCW. However, the mean value of 
k~ is also inversely proportional to the square root of the 
temperature. In regions where the dipolarization is not as 
strong and the temperature is still high the scattering can 
continue. Therefore, the simulation demonstrates that the 
plasma temperature in the SCW can influence when or if the 
proton scattering is quenched during substorm conditions. 

[ 32 ] There are a few features of the substorm that the 
simulation does not reliably reproduce. We have already 
mentioned that the plasma in the near-earth tail was overly 
energetic due to a previous substorm in the simulation. Also, 
the simulation onset latitude occurs right around 70 degrees 
latitude whereas the actual onset latitude was slightly below 
64 degrees. This second deficiency can also be partially 
attributed to the earlier substorm as some of the lobe flux 
was probably transported to the inner magnetosphere via 
reconnection and enhanced convection associated with the 
earlier substorm. That additional flux would inflate the inner 
magnetosphere pushing the transition region to further radial 
distances (i.e., higher latitudes). 

[ 33 ] Another prominent, but less import feature not repro- 
duced by the simulation was the proton aurora extending past 
dusk and well into the afternoon sector. This is most likely 
because the scattering which caused the precipitation was 
most likely due to wave-particle interactions [Fuselier et al., 
2004; Zhang et al., 2008; Spasojevic and Fuselier, 2009] 
which are not included in our model. It could also be the 
result of the magnetospheric compression associated with 
the dynamic pressure increase [Lyons et al., 2005], However, 
the simulation does produce significant electron precipitation 
in that region which is visible above the day-glow in the 
IMAGE-W1C data. 

[ 34 ] Despite the shortcomings of the simulation, it repro- 
duces many of the key features of the substorm (and substorms 
in general) which we now summarize. In the simulation, there 
is an equatorward drifting of the precipitation prior to the ini- 
tial auroral brightening. After the brightening, a westward 
traveling surge develops with the poleward expansion as is 
expected. The simulation produces a second increase in the 
proton precipitation at the same time as a second significant 
increase in the electron precipitation consistent with an 
increase in tail temperature. The observations show that the 
substorm had multiple brightenings (one at 13:07:49 UT and 
one at 13:44:37 UT) and so it is likely that the second bright- 
ening in the simulation corresponds to the second brightening 
in the observations. In the tail, dipolarization is observed in the 
near-earth region driven by bursty bulk flows. Additionally, 
the proton precipitation splits longitudinally in both cases. In 
the simulation and the data, the split region encompasses 
approximately sixty degrees longitude. 


6. Discussion 

[ 35 ] As mentioned briefly in the previous section, from 
equations (1) and (2), it is clear that: 


K 


Rc\B\ 

Vf 


(6) 


Using the nominal CPS temperature and a slightly lower 
\B\ ( T = 4.2 keV, \B\ = 7.5 nT) then the average gyroradius of a 
proton in the growth phase CPS is expected to be approxi- 
mately 0.2 R e . The current sheet thickness is often on the 
order of 2000 km or less [Baumjohann et al, 2007; Zhou 
et al., 2009]. Using these parameters, the growth phase k is 
~ 1 . 2 . 

[36] At the start of the expansion phase, if the field inside 
the SCW dipolarizes completely, it can be shown from 
equation (3) that in the equatorial plane: 

Rc = \ (7) 

where R is the radial distance to the Earth’s center. Thus, 
R c can increase by a factor of ~10 or more at 8 R E . Addi- 
tionally, the magnitude of the magnetic field increases with 
the dipolarization. In the equatorial plane, if the growth 
phase [B| was ~7.5 nT and the field dipolarized com- 
pletely (jR| = |j), the field strength would increase by a factor 
of ^8. These two factors would yield an increase in k by a 
factor of nearly 10 (at 8 R E ) if the field dipolarized com- 
pletely. Of course, n still depends on T~ l/4 . It is common for 
the temperature to increase during a substorm by a factor of 
~4 [Baumjohann, 1991]. Including this temperature increase 
still yields an increase in u by a factor of ~7 at 8 R E . More 
concretely, using equations (2) and (7), the average value for 
k in the equatorial plane of a dipole field is: 

Kdipole = (wvmm) (8) 

where M is the dipole moment (30.4 pT ■ R| for Earth). 
Figure 12 shows equation (8) for a few different plasma 
temperatures. Clearly, in a completely dipolar field, the 
scattering should be quenched. However, if the dipolarization 
is not so complete (if [B| and R c increase by only a factor of 
3.5 for example), but the plasma temperature still increases 
by a factor of 4, k will only increase by a factor of 2.5 which 
results in an average n after dipolarization of ^3 which is 
right at the edge of the scattering threshold (equation (1)). 

[ 37 ] In both of the simulations presented here, the edges of 
the split region in the proton precipitation maps to the 
dipolarized SCW. However, in the April 28, 2001 event, there 
seems to be two levels of dipolarization. The inner, highly 
dipolarized region centered near x = — 6, y = 6 contains the 
majority of the SCW flux. There is also a less dipolarized 
region which extends from about y = 0 to y = 8 (Figures 1 lc 
and llg) which forms a secondary (less dipolarized) SCW 
that encompasses the first. The proton precipitation from the 
outer region is influenced by the high temperature plasma on 
its eastern flank which overcomes the lower curvature to pro- 
duce scattering (Figure 13). 

[38] This suggests that regions of high temperature plasma 
and less stretched but still non-dipolar field can exist during 
substorm conditions. If that energization profile existed 
globally and the field did not dipolarize as significantly, 
there would be no splitting of the proton aurora which is also 
consistent with the speculation by Gilson et al. [201 1], It is 
expected that substorms with that configuration would pro- 
duce weaker field-aligned currents due to smaller pressure 
gradients and smaller flow bursts expected with the weaker 
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Figure 12. Average values for k in the equatorial plane of Earth’s dipole magnetic field as a function of 
distance. For a large range of temperatures, k remains above 3 indicating that it is very difficult to scatter 
protons into the loss cone in a dipolar geometry. 


dipolarization. From this perspective, it is reasonable that the 
proton aurora is less likely to split during weaker substorms. 

[ 39 ] Our model has a few shortcomings that must be 
addressed. The MHD model used in this study cannot 
accurately capture all of the physics associated with sub- 
storm particle injections [e.g., Birn et al., 1997; Spanswick 
et al., 2009]. As these injections tend to accelerate the high 
energy portion of the distribution, they will also make that 
portion of the distribution more likely to scatter and pre- 
cipitate while the lower energy (i.e., thermal) portion of the 
distribution is not available to precipitate. This effect would 
not be picked up by our model and would also decrease the 
likelihood of observing (complete) splitting for a particular 
event. Additionally, magnetotail particle distributions observed 
are typically more like kappa distributions than Maxwellian 
[e.g., Haaland et al., 2010], Since kappa distributions have a 
larger particle flux in the high energy portion of the spectrum 
our calculations probably underestimate the energy flux into 
ionosphere from neglecting those higher energy particles. 

[ 40 ] Even in a thermal distribution of particles, our model 
has some shortcomings. It assumes that the scattering from a 
region of space is determined only by the average gyroradius 
of the distribution. For a given field geometry, this produces 
a critical temperature threshold for turning on and off the 
scattering along that flux tube. In reality, a distribution just 
below the critical temperature and a distribution just above 
the critical temperature will likely produce similar precipi- 
tation. This assumption probably leads to an over estimation 
of the flux into the magnetosphere especially in regions 
where the temperature is close to the critical temperature 
which probably partially offsets the error from assuming a 


Maxwellian particle distribution instead of a Kappa distri- 
bution. In any event, due to these simplifying assumptions, 
the calculated energy flux into the ionosphere is only 
expected to be qualitatively similar to observations. 

[ 41 ] The single fluid MHD used for the simulations pre- 
sented here also does not resolve the physics associated with 
gradient and curvature drifting of ions in the inner magne- 
tosphere. This mechanism is often used to explain a dawn- 
dusk asymmetry in the proton precipitation. However, 
Newell et al. [2009] shows that during moderately strong 
solar wind driving, the asymmetry is not consistent with the 
conventional wisdom. Because of this, it is unclear the role 
that the high energy drifting ions play in the precipitation 
during times of increased geomagnetic activity. In any event, 
those ions are not reproduced in the MHD model and any 
precipitation resulting from those ions is not reproduced in 
our simulations. 

[ 42 ] Finally, the model neglects transit time from the 
magnetosphere to the ionosphere. Since the bounce time for 
a thermal proton is on the order of a minute, this can intro- 
duce some additional timing ambiguity. However, this tim- 
ing ambiguity is generally less than the ambiguity already 
present when trying to interpret the MHD simulation results 
and is therefore probably not significant. 

7. Summary and Conclusion 

[ 43 ] We have presented simulations of the global nightside 
proton precipitation using the output of a global MHD 
simulation and the field line curvature model to determine 
the loss cone filling. The simulation reproduces the global 
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Figure 13. (a-d) Temperature and (e-h) V x on the central plasma sheet. The contours are the same as in 
Figure 1 1. A flow channel opens and drives a dipolarization front into the inner magnetosphere where the 
flux piles up and forms the current wedge. The high temperature region on the eastward flank of the current 
wedge causes the scattering criteria to be satisfied even though the field is marginally dipolarized there. 


features of the observed proton aurora during substorms 
(equatorward motion during the growth phase, intensifica- 
tion and slight poleward expansion at onset) and the general 
ordering of precipitation boundaries presented by Mende 
[2003], The simulation supports the interpretation of 


Sergeev et cil. [1983] that the nightside proton precipitation 
is dominated by field line curvature scattering of CPS pro- 
tons into the loss cone since the simulation captures the 
global features of the proton aurora without including wave- 
particle interactions. While this idea is not new, we have 
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calculated the field geometry and average gyroradius self- 
consistently instead of assuming field geometries and/or 
particle distributions as done previously [Sergeev and 
Tsvganenko, 1982; Sergeev et al., 1983; Delcourt and 
Martin , 1994; Delcourt et al., 1996]. Also, we have com- 
pared our simulation with data from global proton auroral 
images and shown that it reproduces the several key features 
near the substorm onset. 

[ 44 ] Of particular interest to this study is the longitudinal 
splitting reported by Gilson et al. [2011]. The simulation 
supports their interpretation that the splitting is mapped into 
the SOW where the dipolarized magnetic field prevents the 
scattering criteria (equation (2)) from being satisfied. How- 
ever, the simulation for the April 28, 2001 event shows that 
even where the magnetic field is moderately dipolarized, sig- 
nificant temperature (gyroradius) enhancements can still cause 
precipitation. The simulation results suggest that the split 
region always maps into the SCW, but that the edges of the 
split region may not correspond to the edges of the SCW (i.e., 
the SCW may be larger than the region where the split maps). 
In fact, the simulation implies that the edge of the SCW is not 
necessarily an easy boundary to define since there can be a 
highly dipolarized region embedded in a moderately dipolar- 
ized region. This and the energization from substorm particle 
injections are most likely the reasons why not all the substorms 
in the statistical study by Gilson et al. [2011] exhibited clear 
proton auroral splitting. The same reasoning can be used to 
explain why the time delay between substorm onset and pro- 
ton auroral splitting was commonly longer than the proton 
bounce time in the study by Gilson et al. [20 1 1 ] . F inally, in the 
simulation, the location of the IB maps to ~7-8 R E in the 
magnetotail for both substorms examined. 

[ 45 ] In the future studies, we will run the global MHD 
model with a coupled inner magnetosphere model. That 
should improve our mapping and also capture the gradient 
curvature drift effects and more of the particle injection 
physics that the global MHD model lacked in this study. The 
improved mapping and boundary conditions in the inner 
magnetosphere should produce a more reliable estimate for 
the mapping of the IB. Also, since the intensity of the sim- 
ulated proton precipitation depends on the plasma tempera- 
ture, the simulated proton precipitation should enable us to 
develop new tools to use when studying particle injections in 
the coupled global MHD model. Thus, we would like to 
understand the source and location of the peak in proton 
auroral brightness in the magnetosphere. 

[46] Using the inner magnetosphere model of Fok et al. 
[2001], pitch angle distributions will also be able to be sim- 
ulated. With this, temperature anisotropies and growth of 
Electro-Magnetic Ion Cyclotron (EMIC) waves could be 
modeled [Blum et al., 2009]. Using this, it should be possible 
to expand the above described proton aurora code to include 
the precipitation from wave scattering by EMIC waves. 
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